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Abstract 

We discuss how the dijet azimuthal correlations in DIS and real photoproduction at 
HERA probe the differential (unintegrated) gluon distribution in the proton. We find a 
strong dependence of the azimuthal correlation pattern on Bjorken-x, photon virtuality 
and the cut on the jet transverse momenta. A rise of the azimuthal decorrelations is 
observed with decreasing Bjorken-x due to the interplay of perturbative and nonpertur- 
bative effects. We predict a strong rise of the same-side jet rate with photon energy for 
real photoproduction. We discuss conditions for the correlation function to be dominated 
by hard perturbative gluons and ways of constraining the size of the nonperturbative soft 
component. We make some predictions for the THERA energy range. The analysis of 
the energy dependence of the isolated jet and two-jet cross sections in photoproduction 
would be a new way to study the not yet well constrained unintegrated gluon distribu- 
tions and to explore the onset of the pQCD regime. 
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1 Introduction 



Most of the recent studies of \ow-x dynamics relevant for HERA concentrated on the analysis 
of the inclusive structure function, i.e. on the total cross section for 7*p scattering [|TJ and on 
diffraction 0. Because at low x photon-gluon fusion is the dominant underlying mechanism, 
such studies open the possibility of better understanding the so-called unintegrated gluon 
distributions, the quantity first introduced in || . 

The jet studies are known to be a good tool to test perturbative QCD effects. It was 
pointed out already some time ago that dijet production in DIS could be a method to study 
the onset of BFKL dynamics both in photo- M and electroproduction ||. Unfortunately in 
practice, due to unavoidable cuts on transverse momenta of jets, one samples rather large 
values of the gluon longitudinal momentum fraction x g , where it is not completely clear what 
is the underlying dynamics and in particular what unintegrated gluon distribution should 
be used. Up to now mostly real photoproduction data were published by ZEUS and HI 
collaborations 0. Recently the HI collaboration at HERA made available some new results 
for jet production ||. However, in all these analyses only rather inclusive observables have 
been discussed. 

In the present note we discuss the jet production beyond the familiar collinear approxima- 
tion and focus on how more exclusive and more differential jet production observables probe 
the unintegrated gluon distribution. Based on the unintegrated gluon distributions found re- 
cently || from the phenomenological analysis of <j l °* p we explore dijet azimuthal decorrelations. 
In Ref.0 such decorrelation effects were discussed in a perturbative BFKL approach, and as 
a consequence only decorrelation effects at some distance from the back-to-back configuration 
could be analyzed. The extension to the whole azimuthal phase space requires understanding 
the differential gluon densities in non-perturbative soft region. The modelling of the nonper- 
turbative soft gluon exchange as done in Ref . II allows us to extend the region of the azimuthal 
angle between jets closer to tt and to explore the onset of the hard regime. 



2 The formalism 

Here we collect basic formulae used throughout the present note. 

At the parton level the total cross section for quark-antiquark dijet production 7* + p — > 
ji + 32 + A (see Fig.^) can be written in a compact way as: 

*£r id *{z,c?) = fd4>[ a d P i ± f d P i ± Im^i .fi T , L ( X ,Q 2 ,p lt± ,P2 t± ) , (i) 

r l,l,mm r 2,±_,min 

where x and Q 2 are standard kinematical variables. In the formula above f g (x g ,K 2 ) is the 
unintegrated gluon distribution, which will be specified somewhat later, and k is the trans- 
verse momentum of the exchanged gluon. It is related to the quark/antiquark jet transverse 
momenta pi t ± and p2,_L as: 

P2,± = « - Pl,X , = Pi,± + P2,_L + 2pi,±.P2,±COS(j) . (2) 

We have written explicitly lower cuts on the transverse momenta of jets in (TT]). The indices 
T and L refer to transverse and longitudinal photons, respectively. The auxilliary quantities 
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introduced in (|T|) 
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for transverse photons and 



&l{x,Q 2 ,Pi,±,P2,±) 
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Pi, 



(4) 



for longitudinal photons. Above we introduced 



e) = z(l - z)Q 2 + m 2 . 
The unintegrated gluon distribution / 9 is evaluated at 



(5) 



X (. 



M 2 + Q 2 
W 2 + Q 2 



(6) 



where 



Mr 



p 2 ± + m 2 f p\ 
— ! H ~ 



(7) 



is flavour dependent. It is obvious then that at large transverse momenta of jets one samples 
larger values of x g than in the case of total cross section. The scale of the running coupling 
constant in and (f|) is taken to be I 2 = max(K 2 ,£/ +p|)> where for small I 2 the coupling 
constant is frozen as in |J. 

The gluon momentum k is responsible for the jets being not exactly back-to-back in contrast 
to the conventional collinear approximation to leading order. In the following we limit ourselves 
to the region of x 7 ~ 1, where the jets are dominantly produced from the quark box on the very 
top of the gluonic ladder. In doing so we restrict ourselves to leading order parton calculation 
and omit jets from the ladder. 

The principal issue is how the isolated single jet and dijet production samples the uninte- 
grated gluon distribution and when the azimuthal correlation function will be dominated by 
hard perturbative gluons. Although in the present note we address this question based on the 
unintegrated gluon distribution f g (x g , k 2 ) from a recent analysis in || where it was modelled 
phenomenologically to describe structure function data at low Bjorken-rc and the total cross 
section for real 7 + p scattering, we believe that our principal conclusions are to a great extent 
model- in dependent. The following simple two-component Ansatz was adopted in [pf 



f g (x g ,K 2 ) =F soft (n 2 ) 



•FhardipC gi ^ ) 



K 



K 2 + K 2 h 



The parameters k s and Kh determine the scale of the transition from the hard to soft gluon 
region p. 
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The soft nonperturbative component was chosen in the Born form 
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Fso f t(K 2 ) = a soft C F N c ^V{K 2 ) . (9) 

71 [K -f- H so f t ) 



The vertex function V(k, 2 ) is expressed by the single-body isoscalar nucleon form factor 

V(k 2 ) « 1 - F(3k 2 ) . (10) 

The standard dipole parametrization for F was used. 
The hard component was taken in the form 

JW* P , * 2 ) = & B \K 2 f$^ o(Ql - * 2 ) + k 2 M« 2 - Ql) , (ii) 

where is of the Born form (see ||). The not yet specified T v t is calculated from known 
conventional DGLAP parametrizations as 

^ 9 ^ 2 ) = dGDG ^ P{X r K2) ■ (I 2 ) 
a log k z 

The results presented in this note were obtained based on a recent MRST98 LO parametriza- 



tion [12]. As an example in Fig.H we show the resulting unintegrated gluon distribution as a 



function of k 2 for a few different values of x g . The ^-independent soft component is shown 
separately by the dashed line. The model in || is limited to small x g <C 0.1 only, and should 
not be applied for x g > 0.03. For other technical details we refer the reader to M. 

The two-component structure (|j) of the unintegrated gluon distribution as displayed in 
Fig.§ leads to interesting consequences for the dijet azimuthal correlations which will be studied 
in the following section. The hard/soft decomposition of the gluon distribution found in || is 
sufficiently generic to believe that our major conclusions on the onset of the hard regime are 
not strongly model dependent and are relevant also to BFKL driven models. 



3 Results 
3.1 Dijets 

In the present note we shall discuss mainly the effect of azimuthal jet-jet correlations and leave 
the effect of other correlations for a separate analysis. We shall limit ourselves to study the 
region of small Bjorken-x only. The cross section for the dijet production strongly depends on 
cuts imposed on kinematical variables. In order to better demonstrate the effect of coexistence 
of perturbative and nonperturbative effects in the following analysis we shall restrict ourselves 
to cuts on kinematical variables in the so-called hadronic center of mass (HCM) sytem (7*- 
proton center of mass). In the present purposefully simplified analysis we impose the cuts on 
the parton level and avoid extra cuts in the laboratory frame. 

In Fig.^, we present da^p — > j-yjz) jd<\> as a function of HCM azimuthal angle between 
jets for two different values of photon virtuality Q 2 = 4 GeV 2 (left panel) and Q 2 = 16 GeV 2 
(right panel) for a series of Bjorken x. In this calculation, we have restricted the transverse 
momenta of jets to Pi ± M ,P2± M > Pt,cut — 4 GeV and summed over light flavours u, d and s. 
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While in the case of the total cross section the effective mass mj of the quark is responsible 
for confinement effects, in the dijet production the cross section is in practice independent 
of the explicit value of mf. One can observe a strong dependence of the azimuthal angle 
decorrelation pattern on Bjorken x. A closer inspection of both panels simultaneously leads 
to the conclusion that averaging over a broad range of Q 2 would to a large extent destroy 
the effect as it involves automatically averaging over a certain range of x g , the most crucial 
variable for the effect to be observed. A significant part of the effect is due to the interplay 
of the hard(perturbative) and soft(nonperturbative) components. This is explained better in 
Fig. | where the correlation function is decomposed into hard and soft components for two 
rather different values of Bjorken-x: x = 10~ 4 (left panel) and x — 5 • 10~ 3 (right panel). In 
this calculation the virtuality of the photon was fixed to Q 2 = 8 GeV 2 . While at x — 
the hard component dominates, at larger x the soft component becomes equally important. 
In the region of small and small x our results are similar to those from perturbative BFKL 
dynamics ||, which is rather a general feature of the /c^-factorization approach. 

The onset of the hard regime can be best seen in the lower part of Fig.f| where we present 
averaged value of the gluon transverse momentum k (solid line) as sampled at different az- 
imuthal angle between jets. For comparison we also show average transverse momentum of 
the harder (p±,hard = max(px t ±,P2,±), long-dashed line) and softer (p± tSO ft = m i n G°i,_L,P2,±), 
short-dashed line) quark(antiquark) jet. The region of small azimuthal angle (same-side 
jets) samples on average rather large values of gluon transverse momentum k where the per- 
turbative component dominates. Here on average gluon transverse momentum is larger than 
the transverse momentum of each of the two jets. At = 0, < k >=< p± t hard > + < P±,soft >■ 
In our case (x 1 ~ 1) the momentum of the same-side jets k is compensated by the transverse 
momentum of hadrons (minijets) which carry a very small fraction of photon's momentum. 
Close to = 7r the average k is rather small, smaller than the transverse momentum of each 
of the jets. In this region, models for the soft component can be tested. 

In the case of symmetric cut in a sharp (in k 2 ) transition between soft and hard region, to 
a crude approximation, the hard region can be estimated as cos0 > — 1 + n 2 ] /2p 2 cut , where k 
is the transition value in the gluon transverse momentum. In the model in |9[], the border is 
not as well defined and the transition region depends in addition on x g . We wish to mention 
that in general the region of larger (|0 — 7r| ~ 0) is model dependent. In our model, the cross 
section depends on the treatment of the infra-red region i.e. the formula used for a s and the 
choice of its argument. We hope that studying this region experimentally could provide new 
information about genuine nonperturbative effects which are rather poorly known up to now. 

The cut on the transverse momentum is usually a tool to define jets. Studying the depen- 
dence of the result on the cut may however provide some new information. As an example in 
Fig.§ we show the dependence of da /dip on the cut on the HCM transverse momenta of jets 
for fixed value of x = 10~ 3 and Q 2 = 8 GeV 2 and in Fig.|6| the corresponding decomposition 
into soft and hard components for pt, C ut = 2 GeV (left panel) and pt, C ut = 6 GeV (right panel). 
The small pt >C ut = 2 GeV may be slightly academic but is chosen here to better emphasize the 
effect. In both cases we observe the dominance of the hard component at small and the soft 
component near to the back-to-back configuration. In the traditional collinear approximation 
such a separation of soft and hard component is not possible as they are both lumped together 
into the integrated gluon density. In the collinear approximation in first order in a s , jets are 
produced back-to-back and only higher order corrections lead to an azimuthal decorrelation. 
By comparing the two panels in Fig.4 one can observe stronger back-to-back correlations for 
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larger pt lCU t- This is due to the fact that the larger p t region is unavoidably related to larger 
x g , that is to the region of f g in which the soft component is more prominent. 

The experimental identification of the effects discussed here requires good statistics in the 



data sample. Up to now, in practice [IU|, one averages rather over broader range of Bjorken 
x, photon virtuality and jet transverse momenta. Most of the effects discussed in the present 
note are then washed out and the information about the small-x dynamics is to a large extent 
lost. It seems rather difficult at present to study the whole correlation function u>(0) for fixed 
values of x and Q 2 . We suggest that at present, instead of analyzing w(<f>) itself, one could 
study the ratio: 

/ w(4>;x,Q ,p t ,cut) dcj) 

S(x, Q 2 ) = ^ : , (13) 

w(4>;x } Q 2 } pt, cu t) d<f> 



o 

i.e. the percentage of the same-side dijets of all dijets, as a function of Bjorken-x, photon 
virtuality Q 2 and/or the cuts on transverse momenta of jets Pt )(M t- I n Table 1, we show as 
an example the predictions of the two-component model |J for S as a function of Bjorken- a; 
and pt,cut = 4 GeV. Here the virtuality of the photon was fixed at Q 2 = 8 GeV 2 . As seen 
from the table, the two-component model predicts a significant dependence of the same-side 
jet fraction S on Bjorken-x and the value of the lower cut-off on jet transverse momenta in 
HCM. The larger Bjorken x, the lower S and the smaller p^cuti the larger S. In fact these 
two effects are at least partially correlated. It should be noticed that the larger pt, C ut means 
automatically the larger x g (see Eq.(^)). These strong effects predicted here should be easy to 
observe experimentally. However, a more detailed comparison with experimental data should 
include the separation of jets in rapidity/azimuthal angle space [] and/or hadronization effects. 
The higher-order QCD corrections are expected to increase S, especially for larger values of x 
and/or larger cuts on transverse momenta, i.e. in the region where it is small. 

In the model in Ref. || the total (real) photoproduction cross section at energies W < 
100 GeV is dominated by the soft component. Only at very high, not yet available, energies 
the hard component would dominate. At "intermediate" energy available at HERA, the two 
components coexist and their fraction is a smooth function of initial 7p energy. In principle, the 
same stays true for the dijet production and has interesting consequences for the jet azimuthal 
correlations. In Fig.[7|, we show the azimuthal correlation for a few different values of ^p center 
of mass energy W = 50, 100, 200, 500 GeV. We observe strongly rising decorrelation of jets 
when going from fixed target energies to energies relevant at THERA |H| . Similarly as in the 
jet electroproduction, the key variable here is the gluon momentum fraction, which is < x g 
5-10~ 2 at W = 50 GeV and < x g >^ 10~ 3 at the THERA energy of 500 GeV. We get S « 
0.5% at W = 50 GeV. At the THERA energy of 500 GeV we predict about 3.5% of the same 
side jets. The rise of the decorrelation with the photon energy in real photoproduction can 
be traced back to larger values of k sampled by the same-side jets than by the opposite-side 
jets. We note that the rise of the decorrelation with increasing energy obtained here based 
on the model from |J , is rather typical of fc^-factorization approach in general. Therefore the 
experimental confirmation of this effect would be a valuable test of this kind of approach. 



lr rhis is more important for the same-side jets. 
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3.2 Isolated jets in the photon hemisphere 



There is another interesting prediction of the two-component model which we wish to discuss 
briefly now. Up to now we have discussed only cases when both jets have transverse momenta 
larger than a certain cut value pt,cut- I n principle there are also cases (events) with one hard 
ip± > Pt,cut) jet and one soft (p± < Pt, cu t) "jet". In this single jet event x 7 < 1, because 
the transverse momentum of the single quark(antiquark) jet is compensated by a transverse 
momentum of a much softer gluon. In the following we shall call such cases "one jet" events 
for simplicity, as opposed to the previous case called here "two jet" events for brevity. Let us 
compare the rate of such "one jet" events to the previously discussed cases of two hard jets 
in photoproduction. As an example in Fig.|] we compare the cross section for the two cases 
with the lower cut on transverse momentum pt tC ut = 4 GeV. Firstly, we observe that the cross 
section for both cases are of similar order. Furthermore we observe a significantly stronger 
rise of the cross section for the "one jet" case than for the "two jet" case which may be a 
bit surprising at least at first look. This different energy dependence is related to different x g 
and k sampled in both cases. For example for W = 100 GeV and pt, C ut — 4 GeV in the "one 
jet" case < x g >~ 0.01 is substantially lower than in the "two jet" case < x g >~ 0.02. Both 
numbers are, however, substantially larger than average < x g >~ 0.005 sampled in the case 
of total cross section. The effect of k is more complicated as averaged n strongly depends on 
for the "two jet" case. The interplay of the two effects {x g and k) causes the unintegrated 
gluon distributions to be sampled differently in the "one jet" and "two jet" cases. This, 
at least potentially, allows the possibility of a further nontrivial test of unintegrated gluon 
distributions. It would be valuable to compare the present predictions with the predictions of 



standard (collinear) NLO approach (see for instance |T3 



4 Conclusions 

Based on the recent model determination || of the unintegrated gluon distribution in the pro- 
ton we have explored the impact of the soft gluon component and the onset of the perturbative 
regime on the dijet azimuthal correlations. We have predicted a strong dependence of the az- 
imuthal correlation pattern on Bjorken x, photon virtuality and the cut on the jet transverse 
momenta. The effects in the electroproduction could be verified now at HERA, provided a 
careful differential {x ,Q 2 , transverse momentum cut) studies of the dijets are made. 

In real photoproduction we have predicted the rise of the same-side jet rate with the photon 
energy. It would be important to compare the results of the model discussed here with the 
result of the standard collinear approach to understand the potential of such a dijet study to 
shed more light on the low-x dynamics which has been studied up to now in rather inclusive 
processes. Finally, we have found that the study of the energy dependence of the "one jet" 
(defined in the text) cross section would be a new test of unintegrated gluon distributions. 
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Table 1: The fraction of the HCM same-side jets S(x,Q 2 ,p ttCUt ) in % as the function of 
Bjorken-x and the symmetric cut on the jet transverse momenta in GeV for Q 2 = 8 GeV 2 . 
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<t> (deg) <|> (deg) 

Figure 3: The cross section for 7*p — > jij 2 as a function of HCM azimuthal angle between jets 
for Q 2 = 4 GeV 2 (left panel) and Q 2 = 16 GeV 2 (right panel) for several values of Bjorken-x 
= 1(T 4 (solid), 2 • 10" 4 (long-dashed), 5 • 10" 4 (short-dashed), 1(T 3 (solid), 2 • 10~ 3 (long- 
dashed), 5 • 1CT 3 (short-dashed). 
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Figure 4: The decomposition of the HCM azimuthal correlation function w(<f>) into hard (long- 
dashed) and soft (short-dashed) components for x = 1CT 4 (upper left panel) and x = 5-10~ 3 
(upper right panel). In this calculation Q 2 = 8 GeV 2 and p t , cu t = 4 GeV. In the lower panels, 
the average gluon transverse momentum < k > (solid line) and p± tSO ft (short-dashed line) and 
P±Mrd (long-dashed line) are shown correspondingly. 
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Figure 5: The cross section for YP ~ * QQ as a function of HCM azimuthal angle between jets 
for Q 2 = 8 GeV 2 , x = for several values of the HCM symmetric cut on the jet transverse 
momenta: Pt,cut — 2,3,4,5,6 GeV. 
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Figure 6: The decomposition of the HCM azimuthal correlation function w((f>) into hard 
(long-dashed) and soft (short-dashed) components for p t ,cut = 2 GeV (left panel) and p t , cu t = 
6 GeV (right panel). In this calculation Q 2 = 8 GeV 2 and x = 10~ 3 . 
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Figure 7: The correlation function of dijets in real photoproduction for different 7^ center- 
of-mass energies W = 50, 100, 200, 500 GeV. In this calculation a lower cut on both jets 
momenta Pt,cut — 4 GeV was imposed. 
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two jets versus one jet 
P,.„ = 4GeV 




Figure 8: The cross section for jp — > "two jet" (solid) and "one jet" (dashed) cases as a 
function of jp CM energy. In this calculation pt tC ut = 4 GeV. 
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